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I . INTRODUCTION 


The U. S. Standard Atmosphere, 1976 (Ref. 1) has recently 
been published and represents the most current thought on the 
definition of the atmosphere from the surface to 1000 km. The 
representation is that of a steady-state idealized atmosphere 
under moderate solar activity. This report documents a computer 
version of this Standard Atmosphere (US 76). 

The computer version of US 76 was designed to be a 
flexible tool which could be modified to the user's needs with 
minimum effort or knowledge of the mechanics of the program. For 
this purpose the program is segmented into three distinct altitude 
regions which represent separate calculations. The program can 
be used as a single package or may be easily divided into 
pertinent subprograms for specialized use. 

The computer program has been tested and will reproduce 
the U.S. Standard Atmosphere, 1976, table values of each para- 
meter to at least three significant figures. The least accurate 
values arise from calculations requiring numerical integrations, 
and these could be improved by decreasing the tolerance allowed 
in the numerical integrations. 

Comparison of the US 76 with other recent atmospheric 
models reveals a reversal of the trend to attempt to reproduce 
the variability of the atmosphere. It has been noted that the 
structure of the current model lends itself to the same type of 
modifications found in the U. S. Standard Atmosphere Supplements, 
1966 (US 66) which were made to the U. S. Standard Atmosphere, 

1962 (US 62). It would appear reasonable to propose that such 
modification to the computer program presented here would be the 
most natural method for including the variability of the upper 
atmosphere within the context of US 76. 
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Because of the lack of representation of atmospheric 
variability usage of this program is somewhat limited. The 
principal uses of this program are considered to consist of three 
types of user requirements. First, this program provides a 
current and a well documented atmosphere model which could be 
used by researchers requiring representative values for given 
atmospheric parameters in their programs. Second, for comparisons 
with other atmospheric models or atmospheric measurements this 
program could be used for the production of graphical comparisons. 
Third, the program could be used as a substitute for a more 
complex atmospheric model for program development and checkout. 

Much of the code for the US 76 computer program was 
developed and checked out with the use of a programmable hand 
calculator. Consequently, it would be a simple extension of the 
current work to optimize the instruction codes and storage 
requirements of the larger program for use on such calculators. 

One such program is already in use, the calculation of atmospheric 
density from the earth's surface to 32 km has been reduced to a 
56 step program requiring only 35 storage locations on one such 
calculator. With card input available on these calculators it 
would appear that a "portable" U. S. Standard Atmosphere 1976, 
could easily be produced. 

The discussion of the US 76 computer program introduced 
in this section is amplified in the sections which follow. A 
detailed description of the computer program is presented in 
Section II. In Section III the use of the program is described 
along with methods to specialize the program to the user's needs 
without modifying the code. In Section IV the US 76 model is 
compared with other current atmospheric models and some 
suggestions are given as to how the program might be modified to 
account for some atmospheric variability. Finally, the results 
are summarized in Section V. A complete program listing is 
provided in the Appendix. 
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II. PROGRAM DESCRIPTION 


The U. S. Standard Atmosphere, 1976 (US 76) is based on 
an empirical, temperature profile for 45N during moderate solar 
activity (Ref. 1). This temperature profile is shown in Figure 
II-l along with that for the US 62. Below 86 km the pressure 
or density, is obtained from the integration of the hydrostatic 
equation; above 86 km density is obtained from the intergation of 
the diffusion equation for each individual atmospheric constitu- 
ent. The remaining state variable, either pressure or density, 
is obtained from the other two using the gas law for an ideal gas. 

The temperature below 86 km, as shown in Figure II-l, is 
specified by segments which are linear functions of geopotential 
altitude. Between 86 km and 91 km the atmosphere is assumed to 
be isothermal and between 110 km and 120 km the temperature is 
a linear function of altitude. Between 91 km and 110 km the 
temperature is represented by a section of an ellipse. Above 
120 km the temperature increases exponentally and is asyptotic 
to 1000°K. The altitude described is measured along a line 
perpendicular to a gebpotential surface and the altitude is 
taken to be zero along the geopotential surface at mean sea level 
(see Ref. 12, pp 217-219). Reference 1 should be consulted for all 
questions concerning the definition of parameters. 

The US 76 subroutine is divided by temperature regions. 
Figure I 1-2 shows the structure of the subroutine. Subroutine 
US 76 calls one of three subroutines, LDEN, MDEN, or NDEN, depend- 
ing upon the altitude designated by the calling program. These 
three subroutines will calculate the density, pressure and 
temperature within these given altitude ranges using the support- 
ing subroutines and functions shown. The subroutine LDEN has no 
actual supporting routines. The subroutines AUX1 and AUX2 are 
called to calculate other parameters in the 0 to 86 km altitude 
range, such as mean collison frequency and dynamic viscosity. 

Thus, these routines are not called by LDEN but by US 76. 
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Figure I 1-2. Subroutine US 76 












The remaining portion of this section will be devoted 
to descriptions of each routine. The analytical development of 
the equations upon which this program is based are developed in 
Reference 1. Because of the nature of this report it is assumed 
that the reader is familiar with this analytic development or has 
ready access to Reference 1. Thus, the analytic development of 
equations will not be repeated here and equations may be presented 
without a complete description. Input and output parameters 
to these subroutines will be described in Section III. 

2.1 SUBROUTINE US 76 

This routine is the basic calling routine for the U. S. 
Standard Atmosphere, 1976. As such, it determines that the 
altitude requested by the calling program is within the 
appropriate limits; 

0 < Z < 1000 km 

otherwise, all returned values will be zero. Also, this routine 
ascertains what information has been requested by the calling 
program and returns that information (see Section III). 

2.2 SUBROUTINE LDEN 

Subroutine LDEN computes the temperature, density and 
pressure in the region between 0 and 86 km. The temperature is 
computed from the hydrostatic equation and the state equation 
using tabulated values of density for each temperature interval. 
The pressure is calculated from the state equation. 

2.2.1 Subroutine AUX1 

Subroutine AUX1 does not support LDEN but calculates 
the mean air particle speed, the mean free path, the mean 
collision frequency and the speed of sound below 86 km. The 
formulae used in this routine are given in Reference 1. 

2.2.2 Subroutine AUX2 

Subroutine AUX2 does not support LDEN but calculates the 
dynamic and kinematic viscosities as well as the thermal 
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conductivity of the atmosphere below 86 km. The formulae used 
in this routine are given in Reference 1. 


2.3 SUBROUTINE MDEN 


Subroutine MDEN computes the temperature, denisty and 
pressure in the region between 86 and 120 km. In this region 
the temperature is isothermal between 86 and 91 km, an elliptical 
function of altitude between 91 and 110 km, and a linear function 
of altitude between 110 and 120 km. The density is computed 
from the diffusion equation in the form 


T(Z 0> 

n.(Z) = n.(Z 0 ) ex P 



( 1 ) 


where the subscript indicates the atmospheric constituent. The 
term v^yCD^ + K) represents a diffusion transport which is 
parameterized to account for actual constituent distribution 
between 86 and 120 km. This term can be written as 


so that 


V . 
1 

dG(Z) 


D. + K 
1 

dZ 


n.(Z) = 

T(Z ) -G(Z) 

n i (Z 0 ) T(Z) e 

i r z I 

exp \ - / f(Z)dZ/ 

( \ ) 


The integral in Equation (2) cannot be solved analytical- 
ly in this region and must be evaluated numerically. This 
evaluation is performed by the supporting routines INTEG and F(Z). 


The temperature, density , the constituent number density 
and the pressure are all evaluated in MDEN. The density is 
computed from the number densities of the constituents which are 
easily obtained once the integral in Equation (2) is evaluated. 


2.3.1 


Subroutine INTEG 


Subroutine INTEG evaluates the integral 
/* Z T 

U = J f(Z) dZ 
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using Simpson's 1/3 Rule. Because the calculations made in this 
routine are not described in Reference 1, a more complete 
description of this routine will be given in the following 
paragraph. 

Simpson's Rule approximates the function U by 

U k - I < f 0 + 2f l + 4f 2 + 2f 3 + « 4 ♦ . . . + 4f k . 1+ f k ) 

where h is the altitude interval between terms, h = Z^ - 

In programming Simpson's Rule no redundant calculations were 
allowed. Thus, the calculation proceeded in the following 
manner where h = i (Z T - Z Q ) 

I Q = f(Z Q ) + f(Z Q + 2h) 

I ± = f(Z + h) 

1 2 = f(Z Q + h/2 ) + f(Z Q + 3 h/2 ) 

1 3 = f(Z Q + h/4) + f(Z Q + 3 h/4) + f (Z Q + 5 h/4) 

+ f(Z Q + 7 h/4) 

• • • • • 

« • • ♦ • 

• • * • • 

x k - f < z o + 2 FT> + f(z o + jjrr) + • • * 


J i = x o + 4I i 

j 2 = J i - 2Il + 4i 2 

J 3 = J 2 - 2I 2 + 4I 3 


J k = J k-1 " 2I k-l + 4I k 


Then take 


From which the estimates U. are obtained 

k 

= m where: m = h/3 

U 2 = (m/2) J 2 

U 3 = (m/4) J 3 


U k = (m/2 k_1 ) J k 

As programmed the estimate U k is accepted if 


U k - U k-1 

U,_ 


< e 


where e is a predetermined constant, or it is terminated if k 
exceeds some given limit. In the program presented in the 
Appendix 


or 


U k - U k-1 


U, 


< 0.0001 


k 

k = 10 

terminates the numerical integration and U k is returned as the 
final estimate. These values cannot be changed without modifying 
the code. 

2.3.2 Function F(Z) 

The function integrated in INTEG is 


f(Z) = r§ 


g / D i \ r MJC a i # R * , dT~| r3 x 

rItt ’^D i + k ) '[_ m i g dZ J ^ ) 
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which is Eq (36) of Reference 1. In this equation the molecular 
diffusion coefficient D^, the eddy diffusion coefficient K, the 
mean molecular weight M, the acceleration due to gravity g, 
the temperature T, and the temperature gradient dT/dZ are 
functions of altitude. The gas constant R*, the molecular weight 
of the constituent and the thermal diffusion coefficient 
are constants. 

The acceleration of gravity, the eddy diffusion coef- 
ficient, the temperature and temperature gradient are simple 
functions of altitude and can be immediately evaluated from the 
equations given in Reference 1. The mean molecular weight is 
taken to be 28.9644 g/mole between 86 and 100/km; above 100 km 
the mean molecular weight is taken to be 28.0134 for the calcu- 
lation of the number densities of molecular nitrogen, atomic 
oxygen and molecular oxygen. Reference 1 suggests that the 
number densities of Argon and Helium be evaluated above 100 km 
using a mean molecular weight based on the number densities of 
nitrogen and atomic and molecular oxygen at the given altitude. 
Since the mean molecular weight based on the three dominant 
constituents is nearly the same as the actual mean molecular 
weight, the latter was approximated between 100 km and 120 km 
by a polynomial 


m = bQ + b^Z + bgZ 


bgZ' 5 + b 4 Z* 


This value was used only in the calculation of Argon and Helium 
above 100 km. 

The molecular diffusion coefficient is dependent upon 
the total number density. It was impractical to compute this for 
every altitude interval; consequently, the molecular diffusion 
coefficient was approximated by a fourth order polynomial 


D i = c i 0 + C 1,1 Z + c i,2 z2 + c i,3 z3 + c i,4 z4 
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2.4 


SUBROUTINE NDEN 


Subroutine NDEN calculates the temperature, density and 
pressure above 120 km. This routine also calculates the number 
densities of all constituents except hydrogen. 

In the region above 120 km the eddy diffusion coefficient 
K vanishes. Thus, the function f(Z) reduces to 


f(Z) 



T dZ 


(4) 


The functional form of the temperature with altitude allows the 
integral 


U = 


r z 

I f(Z) dZ 


(5) 


to be evaluated analytically. 

2.4.1 Subroutine NHYD 

Above 150 km the hydrogen number density is computed 

in subroutine NHYD. Above 500 km the hydrogen number density 

is computed from Eq (1) with f(Z) defined by Eq (4). However, 

for hydrogen G(Z) = 0 in Eq (2) and the base of the integral 

in Eq (5) is taken as Z = 500 km rather than 120 km. Between 

o 

150 km and 500 km, Eq (1) is modified so that 


n(Z;H) 


n(Z 0 ;H) 



T(Z 0 )\ 
"tTzT ) 


-T 

e 


( 6 ) 
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This is Eq (39) of Reference 1. It should be remembered that 
Z Q = 500 km and 150 £ Z < 500 km when this equation is applied. 
Here 

Z 

x h U = f f(Z) dZ 
Z 0 

The integral in Eq (6) is evaluated numerically by 
Simpson’s 1/3 Rule. The method used is identical to that 
described in Section 2.3.1. Note that if the calculation does 
not converge to the preselected limit, e = 0.0001, the last 
estimate U^q is accepted as the value of the integral and an 
error message is printed. 

2.4.2 Function FX(Z) 

The Function FX(Z ) is given by 



from Eq(6). As before the molecular diffusion coefficient D(H) 
depends upon the total number density. In Function F(Z) the 
diffusion coefficient was approximated by a fourth order poly- 
nomial, here the altitude range for the integral is too large. 
Consequently, the total number density N was approximated by 

tn N = aQ 1 + a 1 i Z + a^Z 2 + a^Z 3 + a^Z 4 

where the superscript is added because the approximation is 
divided into four regions 

1 1 150 < Z < 240 

2 240 < Z < 340 

3 340 < Z < 440 

4 440 < Z < 500 
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The diffusion coefficient was then computed from this total 


number density by 


D(H) = - 
N 



where a, b. and T are constants defined in Reference 1. 

o 
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III. PROGRAM USE 


The US 76 computer program was designed in a modular 
fashion, meaning that certain sections of the program can be 
easily adapted for use without using the entire subroutine 
package. This is accomplished by making various sections 
independent of the other sections of the program and by keeping 
the subroutine interfaces or linkages simple. In this section 
the input and output of the program are described. The common 
blocks are also described so that similar names or parameters can 
be avoided in the user's program or easily changed in the US 76 
subroutine package. The interfaces between routines are 
discussed along with possible segmenting of the subroutine 
package. 

3.1 INPUT/OUTPUT 

Subroutine US 76 has six arguments 

SUBROUT INE US76 ( Z , ID , DEN , PRES , TEMP , A ) 

The only inputs are the altitude Z given in kilometers and an 
indicator ID which is a constant equal to or greater than zero. 

_o 

The outputs are density DEN given in kg • m , the pressure 
PRES given in millibars and the temperature TEMP given in degrees 
Kelvin. Additional output is stored in the vector A which is 
dimensioned A(7) and needs to be so dimensioned in the calling 
program. Below 86 km the vector A may contain dynamic and 
kinematic viscosity, thermal conductivity, mean free path 
collision frequency, particle velocity and speed of sound. All 
of these are given in metric units corresponding to those found 
in Reference 1. Above 86 km the vector A contains the number 
densities of the six constituents and the total number density 
given in per cubic meter. 

The input and output are summarized in Figure III-l. 

The output allows the user to reproduce most of the data found 
in the tables of Reference 1. However, many parameters are not 
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OUTPUT DEN DENSITY IKg/M»« 3I 
PRES PRESSURE <nB) 
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Figure III-l. Input and Output for Subroutine US 76 


made available to the user without suitable calculations or 
conversions made in the calling program. At no altitude is 
geopotential height or any parameter as a function of geopotential 
height given. English units are not used in this subroutine 
package. Ratios of parameters which are found in the tables of 
Reference 1 are not available from these routines. In addition, 
the acceleration of gravity, while calculated in the routines, 
is not available as output nor is the mean molecular weight nor 
is the total number density below 86 km. The pressure scale 
height is not computed. With these exceptions the tables of 
Reference 1 can be reproduced to three significant figures. 

3.1.1 Error Messages 

There are two error messages which the user could 
encounter but are unlikely to occur without some modification 
of the subroutine package or the code itself. In general, the 
only error the user can make, if he calls the US 76 subroutine 
with the proper list, is to specify an altitude outside the 
range. This could occur if a negative altitude or an altitude 
greater than 1000 km. was requested. This latter could easily 
occur if the calling program is working in English units and 
no conversion to metric is made before calling US 76. This error 
will not result in an error message but will result in the 
output for all parameters being zero. 

It should be noted that none of the main subroutines 
LDEN, MDEN or NDEN checks altitude, this is a user responsibility 
if the calling routine US 76 is not used. However, in INTEG 
which is called by MDEN if the altitude is less than 86 km or 
greater than 120 km an error message will be written and the 
integral will be evaluated as zero. The calculation will not 
be terminated. 
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The second error message can be generated either in 
INTEG or in NHYD if the numerical integration fails to converge 
in the specified number of steps. In this case the last computed 
value of the integral is taken and the calculation allowed to 
continue. This contingency has not occurred under normal use and 
would probably only occur should the user modify the convergence 
criteria e 



or modify the number of steps, k, taken before discontinuing the 
numerical integration. 

3 . 2 COMMON BLOCKS 

The only common blocks occur in routines NDEN, INTEG and 
FX(Z) which are called for altitudes above 120 km. Two common 
blocks are defined in NDEN 

COMMON /BLK2/ AH, XN11, Til, Zll , XMH 
COMMON /BLK3/ TINF, RAD, XL. 

These blocks are used to transfer data to NHYD and FX(Z) and 
could be replaced with the appropriate DATA statement. 

3.3 PROGRAM MODIFICATION 

3.3.1 Subroutine Interfaces 

The interfaces between the subroutine US 76 and the 
main subroutines are all similar. The call to LDEN requires 
only four parameters 

CALL LDEN (Z, DEN, TEMP, PRES). 

The altitude Z which is input from the calling routine to LDEN 

and density DEN, temperature TEMP, and pressure PRES which are 

output from LDEN to the calling program. The calling sequence 

for MDEN and NDEN are identical and are the same as LDEN except 

for two additional terms 

CALL MDEN(Z, DEN, TEMP, PRES, A(1),A(7) ) 

CALL NDEN(Z, DEN, TEMP, PRES, A(1),A(7) ) 
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The first of these terms, which as shown here as A(l), is a 
variable, dimensioned at least 6 locations, and is used to out- 
put number densities of the constituents. The other term A(7) 
requires only one location and is used to output total number 
density. 

The routines AUX1 and AUX2 are independent of LDEN but 
require output from LDEN to complete their calculations. Sub- 
routine AUX1 requires 

CALL AUX1(Z, TEMP, PRES, A(6) ,A(4) ,A(5) ,A(7) ) 

altitude, temperature and pressure as input. The output consists 
of particle velocity A(6), mean free path A(4), collison frequency 
A(5), and speed of sound A(7). Subroutine AUX2 requires altitude, 
temperature and 


CALL AUX2 ( Z , TEMP , DEN ,A(1),A(2),A(3) ) 

density as input. The output consists of dynamic A(l) and 
kinematic A(2) viscosity and thermal conductivity A(3). 

The linkage between MDEN and INTEG consists of only 
three parameters. The input to INTEG is the 

CALL INTEG (N,Z ,C) 

constituent index N and the altitude Z. The output is the value 
Of the numerical integration. The interface between INTEG and 
F(Z) is a function 

FUNCTION F(Z,N) 

link requiring as input the altitude and constituent index. 

The link between NDEN and NHYD consists of three input 
parameters and one dimensioned variable to provide output . The 
inputs are the altitude Z, 

CALL NHYD ( Z , T , T10 , XN ) 
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the temperature T and a constant temperature at 120 km T10. The 
output requires a dimensioned variable XN(6) but only the last 
position is modified. The interface between NHYD and FX(Z) is a 
functional linkage requiring only the altitude Z. 

FUNCTION FX(Z) 


as input. 

3.3.2 Modification of Program 

The simpliest modification of the subroutine package 
would be to remove the calling routine as shown in Figure III-2. 
This allows the user to select only that portion of the code 
relevent to his problem. 

If the user is interested only in the 0 to 86 km 
altitude range then AUX1 or AUX2 could be deleted. Alternately 
if the parameters computed by AUX1 and or AUX2 were required for 
known input then they could be used independently. 

Another modification would be the construction of a new 
calling routine to produce the thermosphere model shown on 
Figure III-3. This would be necessary if the entire altitude 
range 86 to above 120 km was required. 

Finally, very simple modifications of the code of 
INTEG would permit the use of this routine as a general Simpson's 
1/3 Rule for which the user supplies the function F(Z). 
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IV. MODEL REVIEW 


The U.S. Standard Atmosphere, 1976, is a static model 
representative of conditions in mid-latitudes during moderate 
solar activity. As such it presents a single temperature profile 
from the surface to 1000 km as well as single profiles of density, 
pressure, mean molecular weight and all other parameters. The 
US 76 is useful for obtaining representative values of atmospheric 
parameters when deviations from these values are not important 
to the user. However, "because of variability of atmospheric 
conditions with spacial location and solar conditions, invariant 
models of the earth's atmosphere (90 to 2500 km) are not useful 
for most engineering applications" (Ref. 2). Moreover, where 
engineering applications tend to require a global view as opposed 
to a mid-latitude view, a similar statement must hold for the 
atmosphere from the surface to 90 km. 

The variability of the atmospheric structure can easily 
be seen from spacial and temporal distributions of the tempera- 
ture. Figure IV-1 illustrates the distribution as observed, in 
the mean, between the surface and 80 km (Ref. 3). Comparison of 
this figure with Figure II-l reveals the similarity of the US 76 
temperature profile with that for the summer months at about 
45° latitude: ~ 290 P K at the surface, 220° at 20km, 280°k at 50 km 
and 190°K above 70 km. A somewhat different profile would be 
indicated even for the winter months at mid-latitudes. Figures 
IV-2 and 3 give a more detailed description of the region between 
40 km and 90 km. From these three figures a general agreement 
about the temperature distribution can be seen even if the 
specifics differ. Finally, Figure IV-4 illustrates the observed 
diurnal and latitudinal variability of the thermosphere tempera- 
ture by means of the global distribution of exospheric 
temperature . 

Since the publication of the US 62 there has been a 
trend toward incorporating the observed variability of the 
atmosphere in Reference Atmospheres, particularly the variability 
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‘latitude 


Figure IV-1. Latitudinal and Seasonal Variation of Temperature with Altitude 
at Solstice (Ref. 3) 




WINTER SUMMER 


Figure IV-2. 


Latitudinal and Seasonal Variation of Temperature 
Altitude at Solstice. (Ref. 4) 


with 
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Figure IV-4. Duinal Variation of Exospheric Temperature as Observed 
by Thomson Scatter (Ref. 6) 


above 90 km. Most notable of these Reference Atmospheres were 
the COSPAR International Reference Atmosphere of 1965 (CIRA 1965) 
and the U. S. Standard Atmosphere Supplements, 1966 (US 66). The 
US 66 included latitudinal and seasonal variability, by means of 
tables, of the atmosphere below 120 km and an analytic represen- 
tation of the atmosphere above 120 km accounting for variations 
due to variable solar activity. More recently, a summarization 
of stratosphere and mesosphere data by Groves (Ref. 5), part of 
which is shown in Figure IV-3, and a revised Thermosphere model 
by Jacchia J71 (Ref. 7) have been combined to form the basis of 
the CIRA 1972. Figure IV-5 illustrates the distribution of exo- 
spheric temperature according to J71. The NASA/MSFC Global 
Reference Atmosphere (Ref. 8) also uses Grove's model of the 
middle atmosphere and a Jacchia Thermosphere Model J70 (Ref. 9). 
However, the Global Reference Atmosphere also includes the Four 
Dimensional World-Wide Atmosphere Model (Ref. 10) to provide 
the atmospheric variability from the surface to 25 km. 

The US 76 viewed from a perspective including recent 
trends in the representation of atmospheric variability appears 
to be somewhat of an anachronism and in need of a new "Supple- 
ment". The thermospheric temperature sturcture of US 76 is very 
similar to that found in US 66 and would lend itself to simple 
modification to incorporate solar variability. Figure IV-6 
suggests how that program modification could be made by incor- 
porating a new subroutine to calculate the exospheric tempera- 
ture. However, this simple approach is also anachronistic in 
that it ignores the trend in thermosphere model development . In 
fact, a more likely approach would be the far more complicated 
integration of Jacchia' s latest thermosphere model J77 (Ref. 11) 
into any new Reference Atmosphere, the thermosphere temperature 
distribution for this model is shown in Figure IV-7 . 

Thus, it would appear that the US 76 while being the 
most recent of the Reference Atmospheres, will find limited use 
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EXOSPHERIC TEMPERATURE DISTRIBUTION AT THE EQUINOXES 
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Figure IV-6. A "Jacchia 1976" Thermosphere Model 








- 90 “ 

b) June solstice. 

Figure IV-7. Duinal Variation of Exosphere Temperature from J77 (Ref. 10) 
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in the engineering community. However, because it is the U. S. 
Standard Atmosphere it will certainly be the standard against 
which all other Reference Atmospheres will be measured. 
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V. SUMMARY 


A computer version of the recently published U. S. 
Standard Atmosphere, 1976, has been developed. The computer 
program has been developed in modular form for easy incorporation 
into the user's program and for easy modification for specialized 
uses. Because of the lack of atmospheric variability the US 76 
will find limited use in aerospace research. However, it will 
become the standard against which all other atmospheric models 
will be measured. 
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(FOR# IS MAINTAIN 

c 

C PROGAM US76 CHECK 

C 


DIMENSION A(7) 

2 • 0 , 

ID « 3 

WRITE (6,9002) 

9002 FOHMAT(?l?) 

DO 100 J» 1 #200 

2 » I ♦ 5 , 

CALL US76(Z,1D#DEN,PRES#TFMF»,A) 

WRITE(6,900V) Z, DEN, PRES, TEMP, (A(I),1»1,7) 

9001 F0RMAT(2X#F6,l#3X,E12 # 4,3X,F7 , .2/9X,7(3X#E12 t 4) ) 
100 CONTINUE 
STOP 
END 

IFOR# IS US76#US76 

SUBROUTINE US76(2#ID#DEN#PHES# TEMP, A) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


US76 IS A CALLING SUBROUTINE FOR THE U.S, STANDARD ATMOSPHERE, 1976 


INPUT Z ALTITUPE(KM) 

ID INDICATOR 

OUTPUT OEN DENSITY (KG/M**3> 

PRES PRESSURE (MB) 

TEMP TEMPERATURE (DEG K) 

A (7) ADDITIONAL VARIABLES DEPENDENT ON ID 


INDICATOR INPUT 


c 

Z,LT,66KM 

ID ■ 

0 

NO ADDITIONAL OUTPUT 

c 


ID » 

1 

VISCOSITY AND THERMAL CONDUCTIVITY 

c 


ID ■ 

2 

MEAN FREE PATH, COLLISION FREQUENCY, MEAN 

c 




PARTICLE VELOCITY, AND SPFED OF SOUND 

c 


ID • 

3 

BOTH THE ABOVE 

c 

c 

Z f GE.86KM 

ID * 

0 

NO ADDITIONAL OUTPUT 

c 

c 


ID « 

1-3 

CONSTITUENT AND TOTAL NUMBER DENSITIES 

c 

c 

INDICATOR OUTPUT 



c 

Z,LT,B6KM 

A(l) 


DYNAMIC VISCOSITY (N*SEC/M**2) 

c 


A (23 


KINEMATIC VISCOSITY (M**2/SEC) 

c 


A (33 


thermal conductivity <J/M*SEC*DEG K) 

c 


A ( 4 ) 


MEAN FREE PATH (M) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


A (5) COLLISION FREQUENCY (1/SEC) 

A (63 PARTICLE VELOCITY (M/SEC) 

A (7) SPEED OF SOUND 

Z.GE.H6KM A(l) N2 NUMBER DENSITY 

A (2) 01 NUMBER DENSITY 

A (3) 02 NUMBER DENSITY 

A (4 ) AR NUMBER DENSITY 

A (5) HE NUMBER DENSITY 

A (6) HI NUMBER DENSITY 

A (7) TOTAL NUMBER DENSITY 
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DIMENSION A (/) 

C 

DEN ■ 0, 

PRES « 0. 

TEMP s 0. 

DO 10 1*1*7 
10 A(I) • 0. 

C 

IF C (Z,LT.0.),ORwCZ.GT.1000,) ) RETURN 

$ 

IF< (Z,GE.0,).AND;CZ.LT,B6 f ) ) 

* CALL LD£N(Z,DLN, TEMP, PRES) 

IE ( (Z,GE.86.).AND.(Z.LT,120.) ) 

* CALL MHENCZ, DEN, TEMP, PRES, A(1),A(7) ) 

IE ( (Z,GE. 120, ). AND. (Z.LE, 1000, ) ) 

* CALL NDEM(Z#DEN, TEMP, PRES, A(1),A(7) ) 

c 

c 

IE ( (I0,EQ,9) ,OR* (ID.GT.i) ) RFTURN 

IFCZ,GE,86,) RETURN 

IE( CID,EQ,t),0Rj(ID.EQ,3> ) 

* CALL AUX2(Z, TEMP, DEN, A(1),A(2),A(3) ) 

IE ( (ID, ED, 2) , OR* (ID,EQ,3) ) 

* CALL AIJXKZ, TEMP, PRES, A (6) , A( 4) , A (5) , A (7) ) 
HFTIJRN 

END 

IFOR* IS LDEN,LPEN 

subroutine LDEN(Z, DEN, T, PRES) 

DATA AE/34, 16319474/ RAD/6356*.7ft6/ 

C 

c COMPUTE geopotential height h 
c 

N * RAD*Z/(RAD ♦ 7) 

c 

c COMPUTE TEMPERATURE 

c 

IE C H t LT. 1 1 , ) 

* T ■ 288,15 - 6 , b * H 

IF c (H f GE!. 1 1 , ) , AND*. (H.LT.20, ) ) 

* T * 216.65 

IF ( (H, GF. 2^.), AND.CH.LT. 32, i ) 

* T * 216,65 ♦ (H - 20,) 

IF t (H # GF. .32, ) , AND*. CH,LT,47, ) ) 

* T « 228,65 ♦ 2;»*(H - 32 f ) 

IE ( <H,Gf„47.) . AND. (H, LT,bl , ) ) 

* T * 2/0,65 

IE ( (H # GF«51.),AND. (H,LT,7i.) ) 

* T « 2/0,65 ** 2,B*( H - 5l t ) 

IEC (H,Gf, 71.), AND.CH.LT, 86,) ) 

* T ■ 214,65 2 . 0 * ( H - 71,) 
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COMPUTE DENSITY 


IF C H.LT.U, ) 

* DEN a I t 225*(288 # l5/T)**n; - AF/6.5) 

IF t (H t GE,n.).AND.(H.LT,2fl,) ) 

* DEN « 3,6392E-01*EXP(-AF*(H - 11,)/216.65) 
IF C (H t GE.20,),ANO.(H.LT,32.) ) 

* DEN ■ 8i803SE~02*(216,65/n**Cr. ♦ AF) 

IF C (H,GE.32,),AND'.(H,LT.47.) ) 

* DEN » 1 ,3225E a *02* (228,65/T) ** ( 1 . ♦ AF/2,8) 
IF t (H,GK,47, ) , AND. (H,l.T,5l , ) ) 

* DEN ■ i.4275E-03*EXP(-AF*(M - 47. ) 7270.65) 
IF C (H,GE,51.) , AND. (H.LT, 71, ) ) 

* DEN a 8,616E»04*(270,65/7)**( 1 , « AF/2',8) 
IF C (H t GE.7l.) f AN0;(M.LT.86,) ) 

* DEN a 6,421 ltf«05* ( 214 ,65/T) ** ( 1’. - AF/2,) 


ORIGINAL PAGF t<s 
OB POOR QUALITY 


COMPUTE PRESSURE 


PRES a 2 .67053072* DEN*T 

RFTURN 

END 


R.fS 4IJX1 . AUXt 

SUBROUTINE AUX1 <Z»T,PRE3,V,Xl, CrCS) 

AUX 1 COMPUTES MEAN AIR PARTICLE SPEED V, MEAN FREE PATH XL 
MEAN COLLISION FREQUENCY C, AND SOUND SPEED CS BELOW 86KM 
GIVEN ALTITUDE Z(KM), TEMPERATURE T(K) , AND PRESSURE PRES(MB) , 

V a 0 , 

XL • 0, 

C ■ 0, 

CS a 0, 

IF C Z,GT.86, ) RETUWN 
COMPUTE MEAN AIR PARTICLE SPEED V(M/8EC) 

V a 27,0.3654*SQRT( T ) 

COMPUTE MEAN FREF PATH XL(M) 

XL • 2,332508F-07*T/PRE3 
COMPUTE MEAN COLLISION FREQUENCY C(1/SEC) 

C a v/XL 

COMPUTE SOUND SPFED CSCM/SEC) 

CS a 20,0468*SURT< T ) 

RETURN 

END 
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[PC.tflS AIJX2,H*Y2 

SUGMIUTXN*' AUX?(2,T,rTN f V7ROVM,VI8KTN#TroNn» 

AUX 2 COMPUTES DYNAMIC VISCOSITY VISDYN, KINEMATIC VISCOSITY VISKIN, 
AND TUP COEPFTCtFNT OF THERMAL CONDUCTIVITY TCOND BELOW 86 KM 
GIVEN ALTITUDE Z(KM), TEMPERATURE T(K) AND DENSITY(KG/M** 3 ) f 

VisnyN « 0, 

VI SKIN * 0, 

TCONO « 0 , 

IF ( Z,GT. 86 . ) RETURN 
COMPUTE DYNAMIC VISCOSITY (N*SEC/M** 2 ) 

VISDYN • t, 458 E- 06 *<T**l # 5 )/(T ♦ 110 , 4 ) 

COMPUTE KINEMATIC VISCOSITY (M** 2 / 3 EC) 

VISKIN « VISDYN/DEN 

COMPUTE COEFFICIENT OF THERMAL CONDUCTIVITY CJ/M*SEC*K) 

TCONO ■ ( 2 , 647 B 2 E- 03 *( T**t ,5 ))/(T ♦ 24 S, 4 * 10 * f **(- 12 ,/T) ) 

C 

RETURN 

ENO 
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|F0H# IS MDEN,MDEN OBlG^^^^ry 

SUBROUTINE M0EN(Z»0EN#T,PHE3,XN#XNT0T) OF POOR QuA^ 1 

~ 'DIMENSION' Q(2#5) > WI2>9) >UT2#9)|XMt5T V 

DIMENSION XN (6) , XN0 (5) 

DATA XN0/l f t30E420,8,6E*16,3,031E*19,l t 351E4l8,7 # 582E*14/ 

DATA 0/0., 0 t ,-S;809644E-04,-3 f 41624RE-»03 f t f 3662t2E»04, 

* 0. » 9,434079E* < 05, 0. ,«*2, 467369E**04 # 0,/ W/0, ,0, ,2,706240£"05 , 

* 5 9 008765E«04, 8',333333E«05, 0, , 8. 333333E-05, 0, ,6,666667E*04, 0,/ 

* U/0 # ,0,,56,90311,97,,86,,e,,06,,0’.,86 t ,0,/ XM/28, 01 34, 15,9994, 

* 31,9988, 39,948,4,0026/ 

OATA GAS/8,31432/ BOLTZ/J.380622E-23/ 
NAMELIST/DCMK/N,/,C,FACT1,FACT2,FACT3,FACT 
XN (6) ■ M, 

XNfOT ■ 0, 

XMM ■ 0, 

COMPUTE TEMPERATURE 
T « 186,8873 

IF ( (Z,Gh«91 . ) , AND, (Z «LT , 1 1 0, ) ) 

IT « 263,1905 - 76.3232* C 1 - ( (Z-91’.)/19 f 9429)**2)**,5 
IF(7,GE t 1 10, ) T*240. ♦ 12*(Z-110.) 

COMPUTE NUMBFR DENSITIES 

DO 100 Is 1 , 5 
N ■ I 

CALL INTF.G(N,Z»C) 

FACT1 * 1 86 , 8673/T 
FACT ■ 1. 

IF ( I ,EQ* 1 ) 00 TO 100 

FACT2 « Q(l,N)/i3'.*W(l,N)3*(EXP(^W(l,N)*(86,-Ua#N))**3T • EXP<«* 

*W(l,N)*(Z-U(lfN))**3) ) 

FACT3 ■ 0.0 

IF ( (Z,Lh,97, ) , AND'. (N.EQ,2) ) FACT3 ■^Q(2,N)/f3 t *W(2,N))*CEXPC- 
*N(2,N)*(U(2,N)"86',)**3) - fcXPC-»W(2,N)*(U(2,N )«Z)**3) ) 

IF C CZ,GT.97,),AND.(N,E0,2) ) 

* FACT3 ■ »1, 1062525 
FACT ■ EXP ( -CFACT2 ♦ FACT3) ) 

100 XN ( I ) s XN0(I)*FACT*FACT1*EXP(-C) 

compute density 

DEN ■ 0, 

DO 200 1*1,5 
XNTOT * XNTDT ♦ XN(I) 

200 XMM s XMM 4 XM(IJ*XN(I) 

XMM « XMM/XNTOT 
PRES ■ XnTOT*ROLTZ*T 
OFN ■ (PRES*XMM)/(GAS*T) 

PRfcS ■ PRES* 1 # E»*fl2 
DEN * DEN*l,E-03 

RETURN 
END 
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(FOR* IS INTEG,INTEG 

SUBROUTINE INTEGCN,Z,U2) 

NAMELIST /ICHK/ N*Z,J#U1,U2 
U2 ■ 0, 

IF C (Z e Lr,86),ORi(Z.GT # 120 # ) ) GO TO 300 
J • 1 

ri » C7*»8().)/2, 

XH • H/3. 

AO ■ F(Z,N) 4 F(86.,N) 

X * 86, * H 
A1 • F(X,N) 

Bl « A0 4 4* At 
Ut ■ XK*H1 
A0 * A 1 

c 

10 JL ■ 2**J 
A l » O, 

c 

DO 100 I»l, JL 
K • on *2 ♦ X 
X ■ 86, + K*H/JL 
100 A 1 » A 1 * F(X»N] 

C 

XK » XK/2. 

d2 * Bt - 2*A0 4 4*A1 

U2 ■ XK*P? 

TEST ■ ABS ( <U2 - U1J/U2 )M00, 

IF(TfcBT # l T,0,010) RETURN 
IF ( J,GT, 1®) GO TO 2O0 
J » J 4 1 
Bl * B2 
A0 « A 1 
U1 ■ l'2 
GO TO 10 
C 

200 NRITE<S»i’001) Ul,U2 
RETURN 

300 WRITE(6»VT02) Z 

RF TURN _ „ „ 

9001 F0RMATCT1INTEGRAL FAILED TO C0NVERGE?/2E1?,4> 

9002 FORMAT (?l ALTITUDE OUTSIDE RANGE7/F12.4) 

END 


IFOR, IS F,F 

FUNCTION F(Z,N) 

DIMENSION C(5,5),B(3),A<5),XMC5) 

DATA RSTAH/8,31432/ B/-7,347J77E402,2.656098E401, 

* «3 f 424545E-0l,l # 94132B5E-GH,-4'.09746hE-06/ 

* C/0,,0,, 0 , ,0,,0 # ,*2.57383bE402 f 1 ,02bH95t40i ,-l',495593E»01 , 
*l,l£)00696F.-03,-2;4956E-06,-2,377458F4l!2 f 1 ,0058961-40 1,*1 , 495594E**0| , 

*r,000696E-03,-2;49b6E-06,-2,701A92r*02,l,0‘>6359E401,-,| # 574121E»0l, 

* 1.0b4529F-03,^2.ft3?74E-06,-2',35900TE4U2;9,253264, 

* -.l t 379937E-01,9 # 279475E«04,-2,327l7(,E^06/ 

DATA XM/28,0134, 1 S', 9994, 31 ,9988, 39', 94 8, 4 ,0026/ A/0,, 0, ,0. ,0,,«.4/ 
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COMPUTE TEMPERATURE SPSP** 14 PAQ E IS 

OF POOR QUALITY 

T ■ 186,8673 

IF((Z t GT.91 # ),ANO;(Z,LT,110,)) 

IT ■ 263,1905 - 76,3232*<1 - ( (Z-9r.)/19 f 9429)**2>**,5 
IF(Z,GE,116,) T i 240, ♦ 1 2* CZ^l 10*3 

COMPUTE GRAVITY 

G ■ 9,80665* (6356*, 766/(6356, 766 ♦ Z) >**2 
FACTi ■ G/(RSTAR*T) 

compute mean molecular weight 

XMM ■ 28,9644 
IF(Z,l.E,100,) GO TO 30 
XMM * 28,0134 
IF(N.GT,3) 

* XMM ■ 8(1) ♦ Z* (B (2) ♦ Z* (B(3) *Z*(B(4) ♦ Z*B(5) ))) 

30 CONTINUE 
FACT2 ■ 1, 

FACT3 a XMM 
FACTA a 0. 

IF(N,EQ«1) GO TO 50 


COMPUTE TEMPERATURE GRADIENT 


DT • 0,0 

IF(<Z # GT,9l,),ANl>;<Z,LT,H0,)) 

* DT ■ -3,8271*((Z<*91,)/C**19,9429))/tl«*C(7-*91 f )/19 # 9429)**2>** # ® 

IF (Z ,GE , 1 10 , ) DT a 12,0 

COMPUTE FDDY DIFFUSION COEFFICIENT 

XK a 0, 

IFCZ.LE.S5,) XK a 120, 

IF C (Z,GT.95.),AN0.(Z,LT,115,) ) 

* XK a 1 20, *EXP ( 1 - (400, /(400, - (Z « 95,)**2)) ) 

COMPUTE MOLFCULAR DIFFUSION COEFFICIENT 

0 ■ EXP ( C ( 1 , N) ♦ Z* (C (2f N) + Z*(C(3,N) * Z*(C(4,N) ♦ Z*C(5,N))))) 
20 CONTINUE 

COMPUTE F(Z) FROM EQUATION 36 

FACT2 a l./(D ♦ XK) 

FACT3 a D*XM(N) ♦ XMM*XK 
FACT4 a (A(N)*RSTAR/G)*DT*D 
50 F a FACTl*FACT2*(FACT3 ♦ FACT4) 

return 

END 
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rP0R#Z8 NDEN, NOEN 

lUBROUTlNE NDEN (Z,T, DENS, PRES, XN, XNTOT) 

U8 76 STANDARD ATMOSPHERE FOR ALTITUDES ABOVE 120KM 

NDEN COMPUTES NUMBER DENSITIES OF ALL NEUTRAL ATMOSPHERIC 

CONSTITUENTS EXCEPT HYDROGEN ABOVE 120KM GIVEN ALTITUDE Z AND TEMPERATURE 


dimension XN (6 J 

COMMON /BUS/ AH, XN11, Til, ZU, XMH 
COMMON /BLK3/ TINF , RAD, XL 


DATA A/0,0,0,0,-e.4/ AH/-tKJ.S5/ T10/360./ T1 1/999,2356/ TINF/1000,/ 
*Z10/120,/ Z1 1/500'./ GAS/H,31 432/ BOLTZ/t .3B06622E-23/ XL/0.01875/ 
*B/0,B60596/ 0/0.0011362/ RAD/«356,7A6/ 

DATA O/0,-5 # t?0»644E-04, 1 .366212E-04, 9.434079E-05, -2 # 457369E*04/ 
*W/0« 2,7fl624>;F*05, 8,333333E-0b, B;3333?3F.*05» 6, 666667E*04/ 

*U/0, 56,9631 1 , BA. , 86., 86,/ 

DATA XN10/3, 726E + 17, 9.275EM6, 4,395E*16, 1,366E*15, 3,888E*13/ 
*XM/28,0134 r 15.9994, 31,9988, 39,946, 4.0026, 1,00797/ XMH/1, 00797 
4/XN11/8.0E4-10/ 


XN (6) « 0. 

COMPUTE GF.OPOTFNTIAL ALTITUDE E EQ(OI) ♦ 2 
E ■ (Zt-Z10) * (RAD * ZIP) / (RAD ♦ Z) 

COMPUTE TEMPERATURE 

T ■ TINF - (TINF - T10)*EXP(-XL*E> 

COMPUTE NUMBER DENSITIES OF ALL CONSTITUENTS EXCEPT HYDROGEN 


DO 50 T«l,5 

TFACT ■ (T10/T)**Ci, ♦ A ( I ) ♦ B*XM(TJ) 

EFACT ■ EXP(-G*XM(I)*E) 

FACT ■ 1,0 
IF(IfFQtl) GO TO 10 
FACT1 ■ EXPC-W(I)*(Z - U(I))**3) 

FACT2 ■ EXP(-W(I)*(Z10 - Ull) )**3) 

FACT ■ EXP( 0(13 * (FACT 1 - FACT2)/(3.*N( 1 ) ) ) 
10 XN(T> * XN10(I)*TFACT*EFACT*FACT 
60 CONTINUE 

COMPUTE TOTAL NUMBER DENSITY 

XNTOT ■ 0. 

DO 200 1*1,5 

200 XNTOT ■ XNTOT ♦ XN(I) 

C 
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COMPUTE HYDROGEN NUMBER DENSITY IF Z GT 150 


ORIGINAL PAGE IS 
OF POOR QUALITY 


IF (Z,GE t 150«) CALL NHYD(Z, T,T10,XN> 

COMPUTE TOTAL NUMBER DENSITY, PRESSURE, MEAN MOLECULAR WEIGHT, AND DENSITY 
PRES ■ 0 

XNTOT ■ XNTOT ♦ XN(6) 

XM0 ■ 0, 

c 

DO 100 I>1, 6 

XM0 • XM0 ♦ XN(I)*XM<I) 

1*0 PRES « PRES ♦ XN(I)*BOLTZ*T 
C 

XM0 • XM0/XNTOT 

DENS ■ (PRES*XM0)/(GAS*T) 

PRES ■ PRES* 1 t E*»02 
DENS ■ DENS*l,E-03 

c 

RETURN 

END 


(FOR# IS NHYD, NHYD 

SUBROUTINE NHYD(Z,T,T1*»XN) 

c 

c NHYD COMPUTES HYDOGCN NUMBER DENSITY ABOVE 150KM 

c 

DIMENSION XN(6) 

* COMMON /BLK2/ AH, XNll, Til, Zll, XMH 

COMMON /BLK3/ TINF, RAD, XL 

c 

DATA B/0 ,0540664/ G/0, 0010137/ 

c 

c 

XN(6) ■ *; 

IF(Z,LT,150,) RETURN 

c 

c 

E • (Z - Z1!)*CRAD ♦ Z11)/(RAD ♦ Z) 

C 

TFACT « rTH/T)**(l, ♦ AH ♦ B* XMH) 

EFACT ■ F.XP(«G*XMH*E) 

XN(6) ■ XN1 l*TFAf.T*EFACT 
IF(Z,GE,500 # ) RETURN 
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CONFUTE HYDROGEN INTEGRAL EQ(39) USING SIMPSON RULE 

U2 a 0, 

J » 1 

H ■ (Z - 500, )/2. 

XK ■ H/3. 

A0 ■ FX(Z) ♦ FX(500,) 

X « 500, ♦ H 
A1 ■ FX(X) 

B1 a A0 ♦ 4,«A1 
U1 • XK*B1 
A0 * Ai 
C 

10 JL ■ 2,**J 
A1 ■ 0, 

C 

00 100 I'l'JL 
k • (i *• n*2 ♦ i 
X ■ 500, ♦ K*H/JL 
100 At ■ A1 + FX(X) 

c 

XK ■ XK/2, 

B2 • R1 •* 2,*A0 +4.*A1 
U2 ■ XK*B2 

TEST » ARS((U2 *»U< )/U2)*100, 

IF (TEST, LT, 0,01 ) GO TO 300 
IF(J,GT,10) GO TO 200 
J ■ J ♦ 1 
B1 a B2 
A0 a At 
U1 * U2 
GO TO 10 
C 

200 WRITE(6#y0Pi ) Z,lM *U2 

300 XN(6) a XN(6) - 1<2*TF ACT*EF ACT 

9001 FORMATC? HYDROGEN CALCULATION DID NOT C0NVERGE74X# ?Z ■ ?F6 # 1,4X,?U 
* 1 *?El? t 4,4X,?U2 a?t!2,4) 

RETURN 

ENO 
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IFOR# IS FX,FX ORIGINAL PAGE1S 

FUNCTION FX(X) OF POOR QUALITY 

COMMON /BLK2/ AH, XNU, Til, Zll, XMH 
COMMON /RLKJ/ TINF, RAD, XL 
DATA 8/0.0540664/ 6/0,0010137/ 

DATA AI/3.305E+21/ BI/0, 8/ Z10/120'. /TO/273,15/ PHI/7.2EM1/ 
DATA T10/360./ 

C 

E » (X - 711)* (RAD + Z1D/CRAD ♦ X) 

EE« (X - 710) * (RAD ♦ Z10)/(RAD * X) 

T ■ TINF - (TINF - T10)*EXP(-XL*EE) 


TFACT ■ (T/Tl 1 )**( 1 , ♦ AH ♦ B*XMH) 

EFACT ■ EXP ( G*XMH*E ) 

IF ( X, IE. 248,) 

* XNT ■ 6, 85832167 37E+01 - X* (4.679689977E-01 «, X* ( 2,790152461E«0 

* 3 - X*(8,032792?9bE-tJ6 - X*8,91 1075174E-09) ) ) 

JF( (X,GT.240,)',AND, (X.Lt.340.) ) 

* XNT ■ 4.94773235E+01 - X* ( 1 ,21 9152248E-01 n X* (4, 1845l8744E»04 

* - X*(7,806663756E-07 - X*b‘.7085b6435E-10) ) ) 

IF C (X.GT.340,)', AND,(X,LE S 440', ) ) 

* XNT ■ 4.520087318E+01 - X*(6.793M62E-02 X* ( l ,626498857E«®4 - 

* X*(2,411b0338E-07 - X*1 .435244756EM0) 1 ) 

IF C X.GT.440, ) 

* XNT » 6,2074742E*01 ~ X* (2. 1 1377105E-0) - X* (6, I85501048E-04 * 

* X* (6, 830390 lbE**87 -X*4,81 193! 816EM0) ) ) 

XNT ■ EXP(XMT) 

C 

DFACT ■ (AI/XNT)*(T/T0)**BI 
C 

FX « TFACT*FFACT*PHI/DFACT 
FX ■ FXM.E + 03 
C 

RETURN 

END 


£U.S. GOVERNMENT PRINTING OFFICE 1978-740-193/114 REGION NO- 4 


45 


